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Abstract 

The universality of intermittency in hydrodynamic turbulence is considered based on a recent model for the 
velocity gradient tensor evolution. Three possible versions of the model are investigated differing in the assumed 
correlation time-scale and forcing strength. Numerical tests show that the same (universal) anomalous relative 
scaling exponents are obtained for the three model variants. It is also found that transverse velocity gradients are 
more intermittent than longitudinal ones, whereas dissipation and enstrophy scale with the same exponents. The 
results are consistent with the universality of intermittency and relative scaling exponents, and suggest that these 
are dictated by the self-stretching terms that are the same in each variant of the model. 

Resume 

Intermittence et universality d'un modele lagrangien des gradients de vitesse en turbulence 3D. Le 

caractere universel du phenomene d'intermittence en turbulence est etudie a partir d'un modele recent regissant 
revolution du tenseur des gradients de vitesse. Trois versions possibles du modele, pour lesquelles les hypotheses 
retenues pour le temps de correlation et l'amplitude du forcage sont differentes, sont analysees. Une integration 
numerique des equations montre que les exposants anormaux des moments relatifs sont les memes pour les trois 
variantes du modele. II est de plus montre que les gradients transversaux de vitesse sont plus intermittents que 
les longitudinaux alors que la dissipation et l'enstrophie se comportent comme des lois de puissance de meme 
exposant. Ces resultats sont coherents avec l'universalite des exposants relatifs et suggerent l'importance du 
terme d'auto-etirement, qui est identique dans les trois variantes du modele. Pour citer cet article : L. Chevillard, 
C. Meneveau, C. R. Mecanique 335 (2007). 
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1. Introduction 



Progress in understanding the small-scale structure of three-dimensional turbulent flow requires the 
study of the velocity gradient tensor Aij = dui /dxj , where u denotes the velocity vector. In incompressible 
flow, A is trace-free, i.e. An — 0. The dynamical evolution of A is obtained by taking the gradient of the 
Navier-Stokes equation: 

dAij = A lk A k | v 8 * Aij (1) 

where d/dt stands for the Lagrangian material derivative, p is the pressure divided by the density of the 
fluid and v is the kinematic viscosity. Neglecting viscous effects and the anisotropic part of the pressure 
Hessian entering in Eq. (1) leads to a closed formulation of the dynamics of the velocity gradient tensor 
known as the Restricted-Euler (RE) equations [1,2]. RE equations predict several phenomena observed 
in various experimental [3,4] and numerical [5] studies of turbulence, such as preferential alignments of 
vorticity and preferential axisymmetric extension. Recently, a system of differential equations describing 
longitudinal and transverse velocity increments has been derived from this approximation and predicts 
non-Gaussian statistics (and in particular skewness) of the components of A [6]. However, in this system 
as well as in the RE equations, the neglect of pressure Hessian and viscous effects leads to singularities 
and precludes the establishment of stationary statistics due to undamped effects of the self-streching 
term. To adress this deficiency of RE dynamics, and based on prior works [7,8,9], a new model has been 
proposed [10]: the Recent Fluid Deformation (RFD) closure. It models both pressure Hessian and viscous 
term entering in Eq. (I). The RFD closure is based on the dynamics and the geometry of the deformation 
experienced by the fluid during its most recent evolution. A de-correlation time scale r entering the 
various closure terms has to be specified, as well as a Gaussian forcing term. It was shown [10] that the 
system with r chosen equal to the Kolmogorov scale and a fixed Gaussian forcing amplitude, reproduces 
stationary statistics with a number of geometric features of the velocity gradient, as well as relative 
scaling exponents of high-order moments, in strikingly close agreement to experimental and numerical 
measurements for real Navier-Stokes turbulence. 

In order to explore the possible universality properties of the model, in this Letter we study two different 
choices for the time-scale and also explore the consequences of varying the amplitude of the forcing term. 
We pose the question of whether the anomalous scaling exponents vary from case to case, and also extend 
the analysis to include scaling properties of the dissipation. 

The model proposed by Ref. [10] begins with a change of variables, expressing the pressure in terms 
of the Lagrangian coordinates X. One may define a mapping Mt ,t between Eulerian and Lagrangian 
coordinates: Mt ,t ■ X e M 3 x e M 3 , where x(X, t) denotes the position at a time t of a fluid 
particle which was at the position x(X, to) = X at the initial time to. The Jacobian matrix of the inverse 
mapping obeys d/dt(dX p /dxi) — —A ki (dX p / dx k )- As argued before [10], for a relative short period of 
time (typically when t — t = t, where r is a characteristic Lagrangian decorrelation time scale for the 
velocity gradient tensor), the solution can be approximated as the matrix exponential of the velocity 
gradient itself, namely (dX p /dxi) = (cxp[— rA]) pi . 

With an Eulerian-Lagrangian change of variables, the pressure Hessian can be written in the following 
way: 

dMt) „, dX p dX q 8 2 p(t) „ Tr(A 2 ) , 

dxidxj ~ dx, dxj dX p dX q ~ Tr^ 1 ) ^ T ,tj ' 1 ' 

In the first approximation, spatial gradients of (dX p /dxi) have been neglected [10]. In the second ap- 
proximation, the short-time solution for (dX p /dxi) mentioned before, the isotropy assumption for the 
Lagrangian pressure Hessian (d 2 p/dX p dX q <~ 6 pq ), and the trace-free condition of A^ have been used. 
Moreover, C T is the short-time Cauchy-Green tensor [10]: C T = e TA e rA . This model (Eq. (2)) can be 
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viewed as a local version of the "tetrad model" [8]. The time-scale r is a model parameter and can be 
viewed as a characteristic time scale of the small dissipative scales of turbulence. Recent material defor- 
mation history can also be used to model other Hessian tensors entering in the Navier-Stokes equation 
(Eq. (1)) such as the viscous part. The resulting Hessian of A is modeled as a friction term and the char- 
acteristic time scale entering in the description is given by the (Reynolds number independent) integral 
time scale T [10]: 

„V»A«-i^>A, (3) 
T 3 w 

which is a stationary version of the model of Jeong and Girimaji [9]. Finally, combining Eqs. (2) and (3) 

into Eq. (1), one obtains a model for the dynamic evolution of A along a Lagrangian trajectory [10], 

1 TviC- 1 ) 



dA 



Tr(A ) 



c: 



A )dt + dW . 



(4) 



T^C" 1 ) T T 3 

A stochastic forcing term c?W has been added to model the combined action of large-scale forcing and 
neighboring eddies. The time evolution of A (Eq. (4)) is thus given by eight independent coupled ordinary 
(or stochastic depending on the forcing dW) differential equations. 

In order to understand the roles played by the pressure Hessian (Eq. (2)) and the viscous term (Eq. 
(3)) in Eq. (4), some analysis can be carried out. Doing so for arbitrary initial conditions A(to) is difficult 
analytically because of the high dimension of the phase space. However, following Ref. [8] one may 
consider the decaying case (with rfW = 0) along a particular direction corresponding to strain with two 
equal positive, and one negative, eigenvalues. Along this direction on the 'Vieillefosse tail', the tensor A 
and the evolution of the relevant eigenvalue A(i) from the model (Eq. (4)) are given by 
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(5) 



The solution of the ODE in Eq. (5) is such that X{t) retains the same sign as A(0). Let us also recall that 
in the RE system, i.e. dA/dt = —A 2 +ITr(A 2 )/3, the time evolution of A is simply given by dX/dt = A 2 
and the finite time divergence is given by the solution A(i) = A(io)/(l — tX(to)) in a finite time l/X(to). 
In our case, we see from Eq. (5) that the anisotropic part of the pressure Hessian acts directly against the 
development of the singularity induced by the self-streching term. Indeed, the coefficient in front of A 2 is 
bounded between [— 1; 1] . Thus this model for pressure Hessian can regularize the finite time divergence 
when A 3> t -1 since then the prefactor of A 2 is close to -1 and the solution of dX/dt = —A 2 tends to zero 
at large times. Further discussions of the regularization along the Vieillefosse tail require specification of 
the time scale r, in particular its dependence on the Reynolds number. 



2. Reynolds number effects, intermittency, and relative anomalous scaling 

In this section, three choices to model Reynolds number changes are considered: (I) constant Kol- 
mogorov time scale (this was the case studied in [10] and for clarity the relevant results will be repeated 
here), (II) local time scale, and (III) variable forcing strength. 

Intermittency is studied by examining the scaling of moments of velocity gradients. As in [10] we 
consider both longitudinal (An, no index summation) and transverse (Ay, i ^ j) gradients. Nelkin [11] 
shows, assuming the relevance of the multifractal formalism in the inertial range, that the relative scaling 
of higher order moments of velocity derivatives should behave as a power law 
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These equations are written for cither longitudinal (i = j, superscript L) or transverse (i j, superscript 
T) gradients. The functions V L (h) and T> T (h) are the longitudinal and transverse singularity spectrum, 
respectively. Imposing (3 = 1 in the inertial range (which is exact for longitudinal velocity increments, 
and a good approximation for transverse ones) leads to ((Aij) 2 ) <~ 1Z e , i.e. finiteness of dissipation is 
recovered [12]. 

To proceed and facilitate interpretation of results, as in Ref. [10] we choose a simple quadratic form for 
the singularity spectrum T> L ' T (h) = 1 — (h — c^' T ) 2 / (2c^ ' T ), with c^' T = | + fcj ' T where the parameter 
' T is called the intermittency exponent (see Ref. [13] for further details). 

Case I - Constant Kolmogorov time scale: First we consider the simplest case in which r is a constant 
value. Because r should scale with the Lagrangian decorrelation time-scale of the velocity gradients, it is 
chosen to be of the order of the Kolmogorov time scale tk [10]. This choice gives an explicit Reynolds 

— 1/2 

number dependence to the model since it that case, tk ~ He -As argued earlier already, from Eq. 
(5), we see that when A 3> r _1 , the model pressure Hessian causes the coefficient appearing in front of 
A 2 to switch from 1 to -1. Thus, it acts to counteract the finite time divergence and causes A to decrease 
in time as 1/t. The viscous part is also very important and can be seen as a very efficient damping 
term with a coefficient which grows exponentially with increasing A's. We have checked numerically that 
for any initial conditions, the system Eq. (5) is such that X(t) — ► when t — > +00. The divergence is 
thus regularized. We have also checked numerically that for any other initial conditions for A, i.e. those 
which cannot be written as in Eq. (5), all components of A evolving under Eq. (4) tend to zero in the 
absence of forcing. Without loss of generality, henceforth all variables will be scaled with the time-scale 
T, i.e. t/T — > t and AijT — > A^. Let us denote by T = t/T the only free parameter. The forcing term 
cfW = G\/2dt is Gaussian and its covariance matrix is assumed to be Reynolds number independent 
(see Ref. [10] for details). The model (Eq. (4)) is solved numerically according to Ref. [10] and stationary 
statistics are obtained. The results are examined from the point of view of intermittency and anomalous 
relative scaling properties. 

To facilitate comparison with the cases considered in II and III, in Fig. l(a-b) we present one of the 
results of Ref. [10]. These results are obtained from numerical integration of the model system over long 
periods of time and the evaluation of moments of various orders. Clearly, intermittency is predicted 
because the K41 line (i.e. of slope p/2, dashed lines)) does not fit the computed results from the model. 
We display (solid lines) various predictions obtained with the help of Eq. (6). The results can be described 
well with the parameters c^ = 0.025 and = 0.040 for the longitudinal and transverse cases, respectively. 
Transverse gradients appear to be more intermittent than longitudinal ones. Also, as shown in Ref. [10] 
longitudinal gradients PDF is skewed. As was stressed in Ref. [10], the intermittency parameters cf 1 and 
c\ are very close to those obtained from experimental data, see Ref. [13,14]. As remarked in Ref. [10], 
however, for values of t/T smaller than about 0.05 (corresponding to a Taylor-based Reynolds number of 
order 300 [15]) the predicted statistics become unrealistic. Still, the fact that a model with only 8 degrees 
of freedom derived directly from the Navier-Stokes equations predicts realistic relative intermittency 
exponents (albeit in a limited range of Reynolds numbers) is quite remarkable. The results raise the 
question of how robust these findings are with respect to other possible choices of the time-scale and 
forcing strengths. This is considered in the next two subsections. 

Case II - Reynolds number dependent local time scale: It has been hypothesized that the dissipative 
scale in turbulence is not constant but fluctuates due to the intermittency phenomenon [16]. Consistent 



with this notion of a local fluctuating cutoff scale, here we choose r(i) = r(7£ e )/y 2Tr(S 2 ) where T is 
a dimensionless parameter, and an unknown function of the Reynolds number. Hence in this case, as 
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Figure 1. Relative scaling properties of higher order moments, i.e. ln(|A;j| p ), as a function of the second order moment 
ln(|j4jj| 2 ). Dashcd-linc represents K41 predictions (i.e. p/2). Solid lines: multifractal predictions (i.e. 3 r (p)/2, see text), 
and o (p = 3), □ (p = 4), A (p=5) and V (p=6). Skewness of longitudinal components is represented with •. Case I 
(a-b): T = TK /T = 0.2, 0.15, 0.1, 0.09, 0.08, 0.07, 0.06 and 0.05 (same value as in Rcf. [10]). Case II (c-d): Lo cal time 
scale r/ v /2Tr(S 2 ) with V = 1.4; 1.35; 1.3; 1.25; 1.2; 1.15; 1.1; 1.05; 1 and 0.95. Case III (e-f): Local time scale l/^/2Tr{S 2 ) 
with increasing variance of the forcing term dW = aGV2dt, with 4<r 2 = 3; 3.5; 4; 4.5; 5; 5.5; 6; 6.5. All parameters listed 
correspond to points going from left to right. 

opposed to Case I where r was held constant, the time-scale fluctuates, being a function the variable 
velocity gradients. When T decreases for a fixed A, the predicted pressure Hessian is closer to isotropy. 
Then the system (4) is dominated more by the quadratic (singularity-inducing) self-stretching term, which 
is what one may expect at higher Reynolds numbers. In this case, for a fluctuating dissipativc time scale, 
the ODE appearing in equation Eq. (4) can be solved exactly and one obtains 

Xnb^ 2 p 4 p 2 p 4 p 

Aft) = r , where a — 4e vT5 — e^m and b = 2e +evT2 >0. (7) 

3aA -et*(-6 2 + 3aAo) 

The long time behavior of the solution depends on the sign of the constant a. For T > -^P In 4 w 0.80, 
A(i) — > at large times. For smaller T's, the solution diverges in a finite time (when A > b 2 /(3a)) and the 
model is unable to regularize the divergence predicted by the self streching term. Therefore, similarly to 
the constant Kolmogorov time-scale option considered in I, this approach appears not to allow reaching 
arbitrarily high Reynolds numbers. The system using different values of T is integrated numerically as in 
I, with the same Gaussian forcing dW — Gy/2dt. Resulting moments of velocity gradients are displayed 
in Fig. l(c-d). As in I, anomalous relative scaling and intermittency is obtained. The solid lines again are 
obtained by using a longitudinal intermittency coefficients = 0.025, the same as the one obtained with 
a constant time scale. The transverse coefficient is = 0.045, also almost the same as that obtained in 
I. 

Case III: Variable forcing strength: Another option to model Reynolds number effects is to vary the 
strength of the stochastic forcing through the term dW = aG\/2dt. The variance a 2 is assumed to 
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Figure 2. Relative scaling properties of higher order moments, i.e. ln{E p ), as a function of the average ln(i?}, where E stands 
for "dissipation" E e or "enstrophy" Eq. The cases (a)-(f) correspond to the same model parameters as in Fig. 1. Various 
orders of moments (E p ) are studied: p = 2 (o), p = 3 (□) and p = 4 (o). Dashcd-linc represents K41 predictions (i.e. slope 
p), solid lines multifractal predictions (i.e. slope p + J r (p), see Ref. [19]). 

increase with increasing Reynolds number based on the notion that the relative strength of forcing (com- 
pared to viscous term) from neighbouring and large eddies increases with Reynolds number. The explicit 
dependence on Reynolds number is unknown but since we use relative scaling, the precise relationship 
with Reynolds number is not needed for the analysis. As a time-scale, in this case we use the local time 

scale, namely r(t) = 1/ \J 2Tr(S 2 ), i.e. the former time scale of Eq. (7) with T = 1 (which is large enough 
to insure regularization along the Vieillcfossc tail). Numerically it is observed that increased a leads to 
increased variance of the velocity gradient components. This is consistent with the view that changing 
the forcing may be viewed as modifying the effective Reynolds number. 

Numerical results for the moments of gradients and quantification of relative scaling are presented 
in Fig. l(e-f) for various strengths a of the forcing. Once again, using the representation of Eq. 6, the 
results agree very well with intermittent exponents consistent with c\ = 0.025 and = 0.040. Thus the 
quantitative predictions of anomalous relative scaling in the model appear to be quite robust with regard 
to details of the regularization time-scale and forcing. 



3. Relative scaling properties of dissipation and enstrophy 

The model (Eq. (4)) can also be used to predict the dynamics and statistics of several norms of the 
velocity gradient tensor, namely E e = 2Tr(S 2 ) (the dissipation divided by viscosity), = 2Tr(fif2 T ) (the 
enstrophy, where i~2 is the antisymetric part of A) and E v = Tr(AA T ) (the "pseudo-dissipation" divided 
by viscosity). It is well known that in turbulent flows, these "dissipation fields" are highly intermittent 
[12,17,18]. In homogeneous and isotropic turbulence, (E t ) — (E^) = (E v ), while higher order moments of 
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these quantities are not similarly linked. We will perform then, in a similar fashion as Fig. 1, a relative 
scaling study of these quantities. 

We present in Fig. 2, similar to Fig. 1, the numerical results for relative scaling of E e and obtained 
from a numerical integration of the model (Eq. (4)), for the three cases I, II and III. Clearly, numerical 
results don't follow K41 predictions (dashed line). It is found that for the all three cases the relative scaling 
properties are the same, i.e. again we observe robustness with respect to how Reynolds number effects 
are modeled. Finally, we observe that dissipation and enstrophy scale the same. The solid line shows the 
multifractal predictions [19,20,21,22], using a unique intermittency parameter fi = 0.25. Similar results 
are obtained when studying the relative scaling properties of the pseudo-dissipation E v , i.e. we obtain 
fim rj 0.25 (data not shown). Thus, in the model the dissipation, enstrophy and "pseudo-dissipation" 
display the same or very similar intermittency exponents. The value of \i s=s 0.25 is in excellent agreement 
with previous numerical and experimental investigations [12,18]. 



4. Summary and Conclusions 

Several options to represent Reynolds number in a Lagrangian velocity gradient model using the newly 
proposed RFD closure [10] have been considered. Numerical integrations of the model show that three 
different ways to represent Reynolds number variations have led to the same (universal) intermittency 
relative scaling exponents. Consistent with data, it has been found that longitudinal intermittency expo- 
nent is of order c% = 0.025 and the transverse one is c% = 0.040. Relative scaling analysis for dissipation, 
enstrophy and pseudo-dissipation shows that all three quantities share the same intermittency coefficient 
(/i w 0.25) in the model predictions. The universality of the results highlights the importance of the self- 
streching term in the dynamics of A. Future work will be devoted to test the robustness of the predicted 
intermittency phenomenon with respect to possible modifications of this term. 

At this stage it may be of interest to recall some prior discussions dealing with the various scaling 
exponents in turbulence. On the one hand, and in agreement with our model results, numerically and 
experimentally the longitudinal and transverse velocity gradient intermittencies were found to be different 
(see for instance Ref. [23] for a recent review on the subject). Conversely (and unlike our model's results), 
they have been predicted to be the same from a field theoretic approach (see Ref. [24] and references 
therein). On the other hand, unlike our model results, enstrophy was found to be more intermittent than 
dissipation in numerical flows [25] (and in experiments, albeit using only a single component of vorticity 
[17]), whereas (in agreement with our model) they are predicted to scale the same at infinite Reynolds 
numbers from simple arguments based on the finiteness of the inertial range pressure spectrum [26], 
or from more systematic irreducible group representations [24]. It is therefore interesting to note that 
the present model provides, for a limited range of Reynolds numbers, an example of dynamics of the 
velocity gradient tensor in which transverse velocity gradients are more intermittent than longitudinal 
ones, whereas dissipation and enstrophy scale the same. Clearly more research is needed to elucidate 
the relationships between various exponents characterizing intermittency in turbulence, and possibly to 
clarify generalizations of the refined similarity hypothesis (e.g. such as in [25]). The present approach of 
using Lagrangian dynamical evolution equations [10] should help shed new light on this long-standing, 
important problem. 
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